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Approximately taking into account the higher-order effects due to multiple-gluon 
initial-state radiation, we extract from the latest Tevatron data of prompt J ftp 
hadroproduction the leading colour-octet matrix elements within the nonrelativistic- 
QCD (NRQCD) factorization formalism proposed by Bodwin, Braaten, and Lepage. 
We find that the matrix elements which describe the formation of J/ ip mesons from 
■ colour-octet cc pairs in the angular-momentum states 2S+1 Lj = x 5o and 3 Pj, with 

J = 0,1,2, which are responsible for the excess of the predicted cross section of 
inelastic J/ip photoproduction over the existing HERA data at high values of the 
inelasticity variable z, are significantly reduced. We conclude that it is premature 
to proclaim a discrepancy between the Tevatron and HERA measurements of in- 
clusive J/tp production in the NRQCD framework. We also consider J/i/j mesons 
originating from the radiative feed down of promptly produced XcJ mesons. 
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1 Introduction 



Since its discovery in 1974, the J/ip meson has provided a useful laboratory for quan- 
titatively testing quantum chromodynamics (QCD) and, in particular, the interplay of 
perturbative and nonperturbative phenomena. Recently, the cross section of inclusive 
J I ijj hadroproduction measured in pp collisions at the Fermilab Tevatron turned out 
to be more than one order of magnitude in excess of what used to be the best theoretical 
prediction [[§, based on the colour- singlet model (CSM). As a solution to this puzzle, 
Bodwin, Braaten, and Lepage [|J proposed the existence of so-called colour-octet pro- 
cesses to fill the gap. The central idea is that cc pairs are produced at short distances 
in colour-octet states and subsequently evolve into physical (colour- singlet) charmonia 
by the nonperturbative emission of soft gluons. The underlying theoretical framework 
is provided by nonrelativistic QCD (NRQCD) endowed with a particular factorization 
hypothesis, which implies a separation of short-distance coefficients, which are amenable 
to perturbative QCD, from long-distance matrix elements, which must be extracted from 
experiment. This formalism involves a double expansion in the strong coupling constant 
a s and the relative velocity v of the bound c quarks, and takes the complete structure of 
the charmonium Fock space into account. 

In order to convincingly establish the phenomenological significance of the colour-octet 
mechanism, it is indispensable to identify it in other kinds of high-energy experiments as 
well. In the meantime, there have been a number of interesting proposals and encouraging 
results concerning possible colour-octet signatures of inclusive J/ip production in e + e~ 
annihilation in the continuum || and on the Z-boson resonance 0, hadronic collisions 
in the fixed-target mode @, B- meson decays f|, and photoproduction in the forward 
direction || and at finite transverse momentum [10|, 11J1 . Although the experimental Z- 



boson-decay data from CERN LEP [|T^] are in agreement with the predicted colour-octet 
signature, they do not exclude the hypothesis that the observed prompt charmonium 
signal is produced by colour-singlet processes alone. 

In the case of inelastic J/ip photoproduction, however, NRQCD with colour-octet 
matrix elements tuned |L3| to fit the Tevatron data |] predicts [0?O] at leading order 
(LO) a distinct rise in cross section as z — ► 1, where z is the fraction of the photon 
energy transferred to the J/ip meson in the proton rest frame, which is not observed by 
the HI and ZEUS |y| collaborations at DESY HERA. This colour-octet charmonium 
anomaly has cast doubts on the validity of the NRQCD factorization hypothesis which 
seems so indispensible to interpret the Tevatron data in a meaningful way. 

This paper is an attempt to rescue the NRQCD approach by approximately taking 
into account dominant higher-order (HO) QCD effects. The basic idea is as follows. The 
predicted excess over the HERA data at z close to unity is chiefly generated by colour- 
octet cc pairs in the states l S Q) 3 Pq, an d 3 i"2 jlOlJlllh where we use the spectroscopic notation 



+1 Lj to indicate the spin S, the orbital angular momentum L, and the total angular 
momentum J. On the other hand, in hadroproduction at the Tevatron, the contributions 
from the colour-octet ^So and 3 Pj states fall off much more strongly with increasing pt 
than the one due to the colour-octet 3 S\ state fL3|, which is greatly suppressed in the 
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quasi-elastic limit of photoproduction [|l(],[n]]. Consequently, the nonperturbative matrix 
elements which are responsible for the colour-octet charmonium crisis are essentially fixed 
by the Tevatron data in the low-pr regime. This is precisely where the LO approximation 



used in Ref. [ 13 1 is expected to become unreliable due to multiple-gluon radiation from 



the initial and final states. In Ref. [|T^], this phenomenon was carefully analyzed in a 
Monte Carlo framework and found to significantly increase the LO cross section. In this 
paper, we perform fits to the latest prompt J/ip data taken by the CDF collaboration 
at the Tevatron incorporating this information [16| on the dominant HO QCD effects. 
We then update the NRQCD predictions for inelastic J/ip photoproduction at HERA to 
find out whether the colour-octet charmonium anomaly persists. We also carry out a LO 
analysis for J/ip mesons originating from the radiative feed down of promptly produced 
Xcj mesons. 

This paper is organized as follows. In Section 0, we specify the theoretical input for 
our analysis. In Section |3], we describe our fitting procedure and present our results for the 
various J/ip and Xcj matrix elements. In Section f|, we present the resulting predictions 
for inelastic J/ip and Xcj photoproduction at HERA. Our conclusions are summarized in 
Section |]. 



2 Theoretical input 

We now describe the theoretical framework for our analysis. We calculate the cross 
section for the inclusive production A + B H + X oi the physical charmonium state H 
in the collision of two hadrons, A and B, in the parton model of QCD, i.e. A and B are 
represented by their parton density functions (PDF's). In the case of photoproduction, 
A is a quasi-real photon, which either directly interacts with the partons inside hadron 
B (direct photoproduction) or fluctuates into a bunch of quarks and gluons, which in 
turn interact with the partons inside hadron B, while the photon remnants give rise 
to hadronic activity in the original photon flight direction (resolved photoproduction). 
Depending on the transverse momentum p? of the H meson, we adopt two different 
pictures to describe its production. For values of the order of the charmonium mass 
Mn ~ 2m c , we consider the formation of a collinear cc pair within the primary hard- 
scattering process (fusion). In the limit pt 3> Mh, only those subprocesses survive where 
the cc pair is created from a single high-energy gluon, charm quark or antiquark which is 
close to its mass shell (fragmentation). It is then useful to describe this two-step process 
as a convolution of the hard-scattering cross section for single-parton production with 
an appropriate fragmentation function (FF). In both cases, the factorization theorem of 
NRQCD Q allows for a systematical treatment of the transition from the cc pair, in 
state n, to the physical H meson. The states n = [ 1, 2S+1 Lj], [8, 2S+1 Lj], where 1 and 8 
indicate the colour multiplicity, span the whole cc Fock space. This formalism comprises 
a separation of short-distance parts, which are amenable to perturbative QCD, from long- 
distance matrix elements iO H [n] \ which must be extracted from experiment. In the case 
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of fusion, the differential cross section of ij — > H + X is decomposed as 

^(ij - H + X) = £ ^ (ij - cc[n] + fc) , (1) 

where the cross sections of the partonic subprocess — > cc[n] + k may be calculated in 
NRQCD as a perturbation expansion in the strong coupling constant since the 



renormalization scale /i is typically set by the transverse mass my = J 4m 2 c + "p\ of the 
H meson. In the case of fragmentation, one has 

D^ H (x, fJ,) = J2 d ^n(x, fi) (0 H [n]) , (2) 

n 

where di^ n (x, /i) gives the probability for the parton i to form a jet that includes a cc pair 
in state n carrying the longitudinal-momentum fraction x. The coefficient di^ n (x, /io) at 
the initial scale /io = 2m c involves only momenta of order m c , and can thus be calculated 
within NRQCD in powers of a s (//o)- The evolution of the FF D^h A*o) up to higher 
fragmentation scales = My is ruled by the timelike Altarelli-Parisi (AP) equations, 
which may be conveniently solved in x space [17] . 



The relative importance of the various terms in Eqs. ([1|) and (0) may be estimated my 
means of NRQCD velocity scaling rules Q. In the limit v — > 0, where v is the average 
velocity of the charm quark in the i/-meson rest frame, each of the nonperturbative 
matrix elements (0 H [n]^ scales with a definite power of v. Thus, Eqs. ([1]) and (|2]) can 
be organized as a double expansion in ot s and v. At leading order in v, Eqs. flU) and (||) 
reduce to the standard factorization formulas of the colour-singlet model, which contain 
just the term referring to the state n of the physical H meson, i.e. [1, 3 5i] and [1, 3 Pj] for 
the J /if) and Xcj mesons, respectively. In the case of the J /if) meson, (O j ^[1, 3 S\\^ is 
related to the nonrelativisitic radial wave function at the origin, Rj/^(0), and may thus be 
extracted from the measured leptonic annihilation rate. The QCD-improved |L8| relation 

iW* - e+n = ( ^[i,3 Sll ) (i _ « , ( 3 ) 

where a is Sommerfeld's fine-structure constant and e c = 2/3 is the fractional charge of 
the charm quark. In a similar way, (0 XcJ [ 1, 3 Pj}) may be determined from the measured 
partial widths of the Xcj- meson decays into light hadrons [TTJ or two photons P0| . The 
leading colour-octet matrix elements of the J /if) meson are \0 J ^[8, ^o]), (O j ^[8, 3 Si]j, 
and (o J ^[ 8, 3 Pj]^, with J = 0,1,2. The leading colour-octet matrix element of the 
Xcj meson is (0 XcJ [8, 3 Si]). Due to heavy-quark spin symmetry, the J-dependent matrix 
elements satisfy the multiplicity relations 



O j ^[8, 3 Pj}) = (2J+ 1) (0 J/ ^[8, 3 iV 
XcJ [l, 3 Pj]) = (2J+ 1) (C XC0 [1, 3 P ] 

XcJ [8, 3 Si\) = (2J+ 1) (C> xc0 [ 8, %} ) , (4) 
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up to terms of relative order v 2 . The scaling with v (and m c ) of the various J/ip and Xcj 
matrix elements is indicated in Tables [I] and 0, respectively. 

The leading colour-octet matrix elements for the J ftp and Xcj mesons have been de- 
termined through fits, based on the fusion |p~3"l , |2Tl and fragmentation |22| mechanisms 
at LO, to CDF data on inclusive charmonium hadropro duct ion [I1J2]. In this paper, we 



perform new fits to the most recent CDF data || taking into account information on the 
HO QCD effects induced by multiple-gluon radiation in a Monte Carlo framework 
and consistently combining the fusion and fragmentation pictures. 

The hadropro duct ion or resolved photoproduction of charmonium at finite px via 
fusion proceeds at LO via the subprocesses gg — > cc + g, gq — > cc + q, and qq — > cc + g, 
where q stands for a light quark or antiquark. The cross sections for the colour-singlet 
states [1, 3 Si] and [1, 3 Pj] may be found in Ref. and those for the colour-octet states 
[8? 1( So]) [8, 3 5i], an d [8, 3 Pj] in Ref. [|l^]. As for direct photoproduction, the contributing 
subprocesses are jg — > cc+g and yq — > cc+q, where the photon can either couple to the q- 
or c-quark line. The cross sections for the colour- singlet and most of the colour-octet states 
are presented in Refs. pi and [ffT|, respectively. The cross section of •yg — > cc[8, 3 Si] + g, 
which is not explicitly listed in Ref. |QJ], is obtained from the one of jg —>■ cc[l, 3 Si] + g 
through multiplication with 15/8. At LO, [!, 3 5i] can only be produced through gg — ► 
cc + g and jg — > cc + g. Furthermore, [ 1, 3 Pj] cannot be generated via 75 — > cc + g and 
*yq — > cc + q, i.e. the direct photoproduction of Xcj mesons is a pure colour-octet process. 

The direct and resolved photoproduction of J/ip and Xcj mesons via fragmentation 
was extensively discussed in Refs. P^j25|p6| . The analysis of resolved photoproduction 
Z5| , |2B| may be readily converted to hadroproduction by replacing the photon PDF's with 



and collected in the 
35 1 i is of 0(a s ), while 



the respective proton PDF's. The relevant partonic cross sections are available at next- 
to-leading order (NLO) in massless QCD [B7J. The di-+ n functions in Eq. (0) pertinent 
to i = g, c, c and n = [ 1, ^i], [1, 3 Pj], [8, 3 S\] were found in Rei 
Appendix of Ref. p5l - Since the gluon is a [8, ^i] state, d g 
£ V- > [8, 1 So] anc i ^g^[8, 3 Pj] require the emission of at least one extra gluon and are thus of 
0{a 2 s ) and beyond. As (e> J /^[8, %}), (p J ^[ 8, ^i]}, and (0 J/ ^[8, 3 P/]) all scale as v\ 
it hence follows that the fragmentation process g — > cc[8, ^i] dominates, while the other 
gluon-initiated colour-octet contributions are marginal. The di-+ n functions for i = c,c 
are of 0(a 2 ), so that the colour-octet contributions due to charm fragmentation are 
suppressed by v 4 relative to the corresponding colour-singlet contributions. In contrast 
to the fusion picture, it is therefore justified to disregard the [8, ^o] and [8, 3 Pj] states in 
the fragmentation picture, as was done in Refs. |2^j25|j26[ . 

We take the renormalization scale \i and the common factorization scale My to be 



M 



f 



where rrix is defined above Eq. @. We define the starting scale /i 



of the FF's as /x = 2m c 



M 



J/i>- 



For our LO analysis, we choose CTEQ4L 



and 



GRV-LO [[30] as the proton and photon PDF's, respectively, and evaluate a s from the 
one-loop formula with = 236 MeV p9| . Whenever we work at NLO, we adopt the 



MS renormalization and factorization scheme and employ CTEQ4M |29 ], GRV-HO |30| , 
and the two- loop formula for a s with A^ = 296 MeV p9| . In the case of fragmentation, 
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we start from the NLO partonic cross sections calculated in the MS scheme with n/ = 4 
massless flavours [27], and modify the factorization scheme for the collinear singularities 
associated with final-state charm quarks so as to match the finite-m c calculation. This 
procedure, which was proposed in Ref. |3l ], is equivalent to the matching between the 
massless- charm calculation in connection with the perturbative FF's of Ref. [0 and the 
massive-charm calculation without FF's. 

Unfortunately, not all ingredients which would be necessary for a fully consistent NLO 
analysis are yet available. In the case of fusion, the NLO corrections to the cross sections 
(da/dt) (ij — > cc + k) in Eq. (|l|) are only known for direct photoproduction in the CSM 
||33|| . Furthermore, in the case of fragmentation, the NLO corrections to the d{^ n functions 
in Eq. (Q) are still unknown. 

In Ref. |]3~3" |, it was found that, under typical HERA conditions, the bulk of the NLO 
corrections to the inclusive cross section of direct J/ip photoproduction in the CSM is 
due to the shift in (O j ^[1*Sx)) which occurs if the QCD correction (K) factor fl8[ 



is 



included in Eq. @, while the genuine NLO corrections arising from the partonic cross 
sections in Eq. (|l|) and the proton PDF's only amount to about 20% in the inelastic regime 
z <; 0.9. Here, z is the inelasticity variable defined as z = Pp-pj/tp/Pp-p^, where p p , p 7 , and 
Pj/^P are the proton, photon, and J /if) four-momenta, respectively. In our recent study of 



inclusive J/ip photoproduction via fragmentation [25|, where we included both the direct 



and resolved channels, we found K factors of about 0.9 and 1.8 for the colour- singlet and 
colour-octet processes, respectively. Guided by the observations made in Ref. we 
do not expect the genuine NLO corrections to the dj_> n functions in Eq. (0), which are 
still missing, to drastically change these findings. Vice versa, we infer from the fragmen- 



tation results of Ref. [25] that, in the case of fusion, the unknown K factors for direct 
photoproduction via the colour-octet mechanism and for resolved photoproduction via 
the colour-singlet and colour-octet mechanisms should also be of perturbative magnitude. 
In the case of direct photoproduction via the colour-octet mechanism, this expectation 
is also supported by a recent NLO calculation of the total cross section ||34|| . In fact, 



at the typical photon-proton CM energy W = 100 GeV, the K factors in the [8, ^o] 
[8, 3 P ], and [8, 3 P 2 ] channels turned out to be as low as 1.26, 1.05, and 1.31, respectively 
[34]| . We conclude that, for inclusive charmonium photoproduction at HERA, the present 
uncertainty in the size of the colour-octet matrix elements is likely to be more significant 
than the partial lack of knowledge of the genuine NLO corrections to the partonic cross 
sections in Eq. ([!]) and the d^ n functions in Eq. (|2|). 

The situation should be very different for inclusive charmonium hadroproduction at the 
Tevatron, especially in the low-pr range, where one expects substantial HO QCD effects 
due to multiple-gluon radiation. Such effects were estimated for the fusion mechanism in 
Ref. [JU by means of the Monte Carlo event generator PYTHIA |55] after implementing 
therein the relevant colour-octet processes, and they were indeed found to be very sizeable. 
The resulting K factor for the px distribution of prompt J/ip mesons may be extracted 
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from Fig. 1 of Ref. |16[ and is conveniently parameterized as 

K M = 2.88 - 2.20 • - l) + 1.04 . - l) 2 . (5) 

It ranges between 4 and 2 for 5 GeV < px < 15 GeV. We assume that this K factor also 
approximately applies to the [1, 3 Si], [8, ^o], [8, 3 Si], and [8, 3 Pj] channels separately. 
Thus, in order to estimate the HO-improved cross section of inclusive J/ip hadropro- 
duction via fusion, we consistently evaluate these channels at LO and in turn multiply 
them by the K factor of Eq. (|5|). On the other hand, our approximate NLO treatment of 
hadroproduction via fragmentation emerges from the corresponding analysis of resolved 
photoproduction [^J2^] by replacing the photon PDF's with proton PDF's. If we extrap- 
olate Eq. (|^) to pt = 20 GeV, we obtain K = 1.72. This result is comparable to the value 
1.45, which we find for the dominant fragmentation channel, [8, 3 Si], using the same value 
for (O j ^[8, 3 Si]^ both at LO and NLO. This indicates that the duality of the fusion and 
fragmentation pictures also carries over to higher orders. 



3 Fit to the data of charmonium hadroproduction 

We now describe our fitting procedure. As experimental input, we use the latest CDF 
data samples of prompt J ftp mesons and of J/ip mesons originating from the radiative 
decays Xcj J/ip + ^ of prompt Xcj mesons 0. The J/ip mesons were detected via 
their J/ip —>■ /i + \i~ decays. The data were collected in pp collisions with centre-of-mass 
(CM) energy ^fs = 1.8 TeV at the Fermilab Tevatron and come as the differential cross 
section da/dpr integrated over the rapidity range \rj\ < 0.6. Each data sample consists of 
11 data points ranging from p T = 5.24 GeV to 18.38 GeV. We adopt the measured values 
of the branching fractions B(xcj J/4> + 7) and B(J/ip — > from Ref. | 36]] . We 



determine {0 J ^[ 1, 3 Si]) from the experimental value of T(J/ip — > £ + £~) via Eq. (|3]). 
Our LO and NLO results are listed in Table |l|. 

We first concentrate on prompt J/ip hadroproduction. In Refs. fT3| , |2~l"|| , the J/ip colour- 
octet matrix elements were fitted on the basis of the fusion picture at LO, whereby the 
[8, 3 S\] component was supplemented with leading logarithms from g —>■ cc[8, 3 Si] frag- 
mentation. This was achieved in an ad-hoc way, namely by multiplication with the ratio of 
the corresponding gluon-fragmentation cross sections with and without LO AP evolution. 
Here, we adopt a more rigorous approach. We introduce a variable separation cut p™*, 
which we choose to coincide with any of the 11 CDF data points, as a demarcation between 
the fusion and fragmentation pictures. In the first step, we determine (O j ^[8, 3 S\\/ by 
fitting the data points for pt > p™* in the fragmentation picture. This result is then used 
as input for the second step, in which (o J ^[ 8, l S$\/ and (^O j ^[8, 3 Po]) are determined 
by fitting the data points for px < p™* in the fusion picture. As was noticed in Ref. |13| 



the cross sections of cc[8, ^o] and cc[8, 3 P ] production have very similar p T dependences, 
so that treating (^O j ^[8, ^o]^ and (o J ^[ 8, 3 Po]) as independent fit parameters would 
yield highly correlated results. It was therefore suggested to fix a specific linear 
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Table 1: Values of the J ftp matrix elements and of r resulting from the LO and HO- 
improved fits to the CDF data §. is defined in Eq. ©. The first 10 (last 2) data 

points are described in the fusion (fragmentation) picture. 





LO 


HO 


scaling 




(7.63 ± 0.54) • 10- 1 GeV 3 


(1.30 ±0.09) GeV 3 


[m 3 f 3 ] 




(3.94 ±0.63) ■ 10~ 3 GeV 3 


(2.73 ±0.45) ■ 10~ 3 GeV 3 


[m 3 f 7 ] 




(6.52 ± 0.67) ■ 10~ 2 GeV 3 


(5.72 ± 1.84) ■ 10~ 3 GeV 3 


[m 3 f 7 ] 


r 


3.47 


3.54 




Xdf fus - 


5.97/10 


2.23/10 




Xdf fra - 


1.53/2 


1.73/2 




Xdf tot - 


7.49/12 


3.96/12 





combination of (o J ^[ 8, ^o]) and (O j ^[8, 3 Po]) instead. Following this suggestion, we 
define 

M}'* = {O j I*[% %]) + ^ (<D J '*[8, 3 P ]) , (6) 

where r is to be chosen in such a way that the superposition of these two channels is 
insensitive to precisely how they are weighted relative to each other. We determine r to 
be the ratio of the fit result for (o J W[&, %)) under the condition (p J ^[ 8, 3 P ]) = 

to the one for (p J ^[8, 3 P ]) /m 2 under the condition (O j ^[8,%}) = 0. We fit by 
analytically minimizing the x 2 values. We measure the overall quality of the fit by adding 
the two x 2 values achieved in the determination of (O j ^[8, 3 Si]\ in the upper p T range 

and of in the lower range. Finally, we determine pjP by minimizing the total 

X 2 ■ It turns out that the best results are obtained if the first 10 (last 2) data points 
are described according to the fusion (fragmentation) picture, i.e. the transition between 
these two pictures happens somewhere around pt = 15 GeV. 

Figures |l|(a) and (b) illustrate for the LO and HO-improved analyses, respectively, 
how the fusion and fragmentation cross section compare with the experimental data and 
how they are decomposed into their colour-singlet and colour-octet components. The 
fragmentation results are only shown for p? > 10.91 GeV, while the fusion results are 
displayed over the full pr range. In each case, the total result (solid lines) is obtained 
as the superposition of the colour- singlet (dotted lines) and colour-octet contributions. 
The [8, 3 Si] contributions (dashed lines) of the fusion and fragmentation pictures almost 
coincide. The colour-octet contribution proportional to M//^ (dot-dashed line) is only 
included in the fusion picture, while it is neglected in the fragmentation picture for rea- 
sons explained in Section ||. Notice that the LO [1, 3 Si] contributions in the fusion and 
fragmentation pictures arise from very different sources. On the one hand, the LO fu- 
sion subprocess gg — ► cc + g has no counterpart in the fragmentation picture. On the 
other hand, the LO fragmentation processes g — >• cc + gg, c — > cc + c, and c ^ cc + c 
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correspond to HO contributions in the fusion picture. This explains why the dotted lines 
in Figs. |l|(a) and (b) do not match. As is by now common knowledge, the colour-octet 
processes are indeed necessary in order to reconcile theory with experiment as far as in- 
clusive hadroproduction of prompt charmonium at the Tevatron is concerned. The results 
and the various x 2 values of the LO and HO-improved fits are summarized in Table [1|. 
The errors quoted in Table [I] only include the experimental errors on the CDF data. 
In addition, there are the usual theoretical uncertainties due to the dependence on the 
renormalization and factorization scales and on the proton PDF's, which have been care- 
fully estimated in Ref. [21). Our LO results should be compared with those obtained in 
Refs. [[U|[n|]. We find reasonable agreement for (and r), which was determined to 

be = (6.6 ± 1.5) • KT 2 GeV 3 in Ref. @ and = (4.38 ± 1.15) • 1(T 2 GeV 3 

(for CTEQ4L |2|) in Ref. |TJ. On the other hand, our result for (O j ^[8, 3 Si]) is some- 



what smaller than the value (6.6 ± 2.1) • 10 3 GeV 3 of Ref. fl3| and considerably smaller 



than the value (1.06 ± 0.14) • 10" 2 GeV 3 of Ref. PJ. We attribute this deviation to 



the different implementations of the fragmentation contribution. As we pass from LO 
to our approximate HO implementation, we observe that (O j ^[8, 3 Sx\\ is only slightly 

decreased, by 31%, while M^'^ drops off by more than one order of magnitude. The 
change in (O j ^[8, 3 5i]^> precisely compensates the enhancement due to the K factor in 

the g — > cc[8, 3 Si] fragmentation channel at high pt- The dramatic reduction of M^^ is 
due to the interplay of two effects. On the one hand, the [8, 3 <Si] fusion channel, whose 
matrix element is fixed by the fit in the fragmentation regime, is strongly enhanced at 
low pt due to multiple-gluon radiation according to Eq. @ and only leaves little room for 
possible contributions due the [8, 1 Sq] and [8, 3 Pj] fusion channels. On the other hand, the 
latter channels are also strongly enhanced at low pt by the K factor of Eq. (|5|), so that 
is even further reduced. The situation is nicely illustrated if we compare the relative 
importance of the various colour-octet fusion contributions in Figs. |l](a) and (b). At LO, 
the combined [8, ^o] and [8, 3 Pj] contributions exceed the one due to the [8, 3 Si] channel 
way up to pt = 10 GeV, while, at HO, the [8, 3 Si] contribution is dominant over the full pt 
range considered. Our HO-improved results for (O j ^[8, 3 S\]^ and M; J ^ are very similar 
to those extracted in Ref. |L6| from the previous CDF data with the 1994 GRV-HO 
proton PDF's |37| , which correctly account for the low-x behaviour of the proton structure 
function F 2 (x,Q 2 ) measured at HERA, namely iO J ^[ 8, 3 Si]) = (3.4 ± 0.4) • 10~ 3 GeV 3 

and = (6.0 ± 1.2) ■ 10~ 3 GeV 3 . Finally, we note that the overall x 2 P er degree of 

freedom (xdf) * s significantly reduced, from 0.62 to 0.33, as we pass from LO to HO. 

In the remainder of this section, we discuss the extraction of the Xcj matrix elements 
through a fit to the CDF data on prompt Xcj mesons 0. For simplicity, we work at LO 
in the fusion picture. The fit results for (O Xc0 [ 1, 3 P ]) and (C Xc0 [8, 3 Si}) are displayed in 
Table|2|. The quality of the fit is rather high, with Xbf = 0.35. Our value for (O Xc0 [ 1, 3 P ]) 
is somewhat larger than those extracted from the partial widths of the Xcj- m eson decays to 
light hadrons and two photons P0fl . In turn, our value for (C Xc0 [8, 3 5*i]) is smaller than 
the one found in Ref. The interplay of the colour-singlet and color-octet processes 
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Table 2: Values of the Xcj matrix elements resulting from the LO fit to the CDF data 0. 
All 11 data points are described in the fusion picture. 





LO 


scaling 


(O**[l*P ]) 
Xdf fus - 


(2.29 ± 0.25) ■ 10- 1 GeV 5 
(6.81 ± 1.75) ■ 10- 4 GeV 3 
3.82/11 


[m^v 5 ] 
[m^v 5 ] 



and the comparison of their superposition with the CDF data is illustrated in Fig. [|. 

4 Predictions for charmonium photoproduction 

Having extracted the leading J /if) and Xcj colour-octet matrix elements from the latest 
CDF data taking into account information on HO QCD effects, we are now in a position 
to explore the phenomenological consequences for inclusive charmonium photoproduction 
in ep collisions at HERA. HERA is presently operated in such a way that E e = 27.5 GeV 
positrons collide with E p = 820 GeV protons in the laboratory frame, so that approxi- 
mately \/s — 300 GeV is available in the CM system. In the case of photoproduction, 
the positrons act as a source of energetic, quasi-real photons, whose energy distribution 
is well described in the Weizsacker- Williams approximation by a well-known formula, 



which may be found e.g. in Eq. (17) of Ref. p5 



The HI [Q and ZEUS Jl5| collaborations recently presented their measurements of 
inelastic J /if) photoproduction. In both experiments, the scattered positron was not 
tagged, so that the maximum photon virtuality was Q^ax = 4 GeV 2 . The measured ep 
cross sections were converted to averaged 7p cross sections by dividing out the photon- 
flux factor, which was evaluated by integrating the Weizsacker- Williams formula over 
the considered photon-energy range. The experimental information was presented in the 
form of distributions differential in the photon-proton CM energy W, the J/if> transverse 
momentum pt, and the inelasticity variable z. We adhere to the kinematic cuts used 
in the ZEUS analysis, namely 0.4 < z < 0.8 and p\ > 1 GeV 2 for the W distribution, 
0.4 < z < 0.9 and 50 GeV < W < 180 GeV for the p\ distribution, and p\ > 1 GeV 2 and 
50 GeV < W < 180 GeV for the z distribution. The photon-flux factor corresponding to 
this W interval is 0.106. By contrast, the HI data refer the W interval 30 GeV < W < 
150 GeV, for which the flux factor is 0.150. Actually, the ZEUS measurement only reached 
down to z > 0.4 and was extrapolated to z = by Monte Carlo simulations yielding an 
enhancement factor of 1.10 fl5| . The HI data were treated in a similar way. We undo 



this artificial extrapolation and divide the published W and p\ distributions [0,0 by 
this factor. The contribution due to if)' mesons with subsequent decay into J /if) mesons 
was not subtracted from the ZEUS data. As in Ref. |15| , we thus multiply the theoretical 



predictions by an overall factor of 1.15 to approximately account for this contribution. 
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The data are mostly concentrated in the \ow-px range, where the fusion picture should 
be valid. 

In Figs. |3](a)-(c), we compare our LO predictions for the W, p?, and z distributions, 
respectively, with the ZEUS and HI data. The circumstance that (O j ^[&, l So\j and 

(O j ^[&, 3 -Po]) are not separately fixed by the fit to the CDF data induces some uncertainty 
in the colour-octet contributions to the cross sections of direct and resolved photoproduc- 
tion and thus also in the total cross section. This uncertainty is encompassed by the results 
for (O j /*[8,%]) = and 8, 3 P ]) = and those for (O j ^[ 8, %}) = and 

(p J W[&, 3 P ]) = {m\lr)M J r ^, which are actually shown in Figs. |(a)-(c). From Fig. |(c), 
we observe that the colour-octet contribution of direct phot opro duct ion is dominant for 
z <;0.5. Thus, it also makes up the bulk of the cross sections shown in Fig. ||(a) and (b), 
which are integrated over 0.4 < z < 0.8 and 0.4 < z < 0.9, respectively. As is evident 
from Figs. [|(b) and (c), this contribution is also responsible for the significant excess of 
the LO predictions over the experimental results for dojdp\ and da/dz at lowpr and high 
z, respectively. (By contrast, such an excess does not show up in Fig. ||(a), where the z in- 
tegration only extends up to 0.8.) This observation was assessed in the ZEUS publication 
|i~5| by the statement that the predictions of a specific leading- order colour-octet model, 
as formulated to describe the CDF data on J / ip hadroproduction, are not consistent with 
the data. To our mind, the salient point is that those predictions were made on the basis 
of LO fits to the CDF data. It is crucial to notice that the colour-octet contribution of 
direct photoproduction is essentially due the [8, ^So], [8, 3 Po], and [8, ^P-^ channels [p~0| , p~lT| 
and thus proportional to M/'^. In Section we have seen that the fit result for is 
greatly reduced by the inclusion of HO QCD effects on hadroproduction at low px Jl6| . 
This nourishes the hope that the apparent discrepancy between theory and experiment 
at HERA may thus be alleviated at higher orders. 

In Figs. El(a)-(c), we take a first step towards a full NLO description of the HERA data 
on inelastic J/ip photoproduction by repeating the analyses of Figs. |3](a)-(c) using the 
HO-improved results of Table [l], the NLO proton and photon PDF's specified in Section |2|, 
and the two- loop formula for a s . However, we do not include the NLO corrections to the 
partonic cross sections in Eq. (HD As already mentioned in Section ^, they are unknown, 



except for direct photoproduction in the CSM [53]. In the latter case, they were found 



to moderately increase the cross section, by about 20% p3| . We observe that the HO 



improved predictions tend to undershoot the data leaving room for a substantial K factor 
due to the missing NLO corrections to the partonic cross sections in Eq. (Jl]). Now, 
the colour-singlet contribution of direct photoproduction, which is well under theoretical 



control [|33j, is by far dominant, except in the corners of phase space, at £<^0.15 and 
z ^ 0.85, where the colour-octet contributions of resolved and direct photoproduction, 
respectively, take over. It has been noticed P3|,EB|,P3 that special care must be exercised 



in the treatment of the latter contributions. In the vicinity of the upper endpoint, HO 
nonperturbative contributions are enhanced and lead to a breakdown of the NRQCD 
expansion, so that only an average cross section over a sufficiently large region close to 
z = 1 can be predicted P9| . On the other hand, one my infer from the study of resolved 
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J /if) photoproduction via fragmentation ]^5],^6| that the corresponding partonic cross 



sections in Eq. (jl]) are likely to receive significant NLO corrections at low z. Of course, 
we should also bear in mind that the predictions shown in Figs. |3| and ^ still suffer from 
considerable theoretical uncertainties related to the choice of the scales fi and Mf, the 
PDF's, and other input parameters such as m c and A*- 4 -* pl| , |33|j . From these observations, 
we conclude that it is premature at this point to speak about a discrepancy between the 
Tevatron and HERA |I^Jl5[| data of inclusive J / if) production within the framework 



of NRQCD §. 

Finally, we present LO predictions for the inclusive photoproduction of Xcj mesons 
with subsequent radiative feed down to J /if) mesons. We work at LO in the fusion picture 
using the Xcj matrix elements of Table [2| In Figs. |5|(a)-(c), we present our results for the 
W, pip, and z distributions, respectively, under the same kinematical conditions as in the 
J / if) case considered in Figs. ^| and f|[ Notice that direct photoproduction of Xcj mesons 
is prohibited at LO in the CSM. From Fig. |§(c), we observe that, except for z^0.8, 
the dominant contribution arises from resolved photoproduction in the CSM. Comparing 
Figs. ^|(a) and (b) with Figs. |3|(a) and (b), we see that, after integration over z, this 
indirect source of J / if) mesons is suppressed by roughly two orders of magnitude relative 
to prompt J /if) production. 

5 Conclusions 

We determined the J /if) colour-octet matrix elements which appear in the NRQCD expan- 
sion Q at leading order in v by fitting the latest Tevatron data of prompt J / if) hadropro- 
duction @. We found that the result for the linear combination of ( O J ^ [ 8, 1 So] ) 

and \O J ^[8, 3 P ]) is substantially reduced if the HO QCD effects due to the multiple 



emission of gluons, which had been estimated by Monte Carlo techniques [16], are taken 
into account. As an important consequence, the intriguing excess of the LO NRQCD 



prediction for inelastic J /if) photoproduction at z close to unity [|K],|ll|] over the HERA 



measurements disappears. We assess this finding as an indication that it is pre- 



mature to proclaim an experimental falsification of the NRQCD framework on the basis 
of the HERA data. Although we believe that our analysis captures the main trend of 
the HO improvement, we stress that it is still at an exploratory level, since a number of 
ingredients which would be necessary for a fully consistent NLO treatment of inclusive 
J /if) hadroproduction and photoproduction are still missing. 
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FIGURE CAPTIONS 



Figure 1: (a) LO and (b) HO-improved fits to the CDF data on the inclusive hadroproduc- 
tion of prompt J/ip mesons [||, which come in the form of da/dpr integrated over \rj\ < 0.6 
as a function of pt- The fragmentation results are only shown for pt > 10.91 GeV, while 
the fusion results are displayed over the full px range. 

Figure 2: LO fit in the fusion picture to the CDF data on the inclusive hadroproduction 
of prompt Xcj mesons with radiative feed down to J/ip mesons ||, which come in the 
form of da/dpr integrated over \r)\ < 0.6 as a function of pt- 

Figure 3: The ZEUS data on the inclusive photoproduction of prompt J/ip mesons fT5| , 
which come in the form of (a) a integrated over 0.4 < z < 0.8 and p\ > 1 GeV 2 as a 
function of W, (b) da/dp? integrated over 0.4 < z < 0.9 and 50 GeV < W < 180 GeV 
as a function of p\, and (c) da/dz integrated over p\ > 1 GeV 2 and 50 GeV < W < 
180 GeV as a function of z, are compared with the respective LO predictions. The total 
results (solid lines) are built up by the direct-photon colour-singlet (dotted lines), direct- 
photon colour-octet (dashed lines), resolved-photon colour-singlet (dot-dashed lines), and 
resolved-photon colour-octet (dot-dot-dashed lines) contributions. The respective HI data 



14| , which refer to 30 GeV < W < 150 GeV, are also shown. 
Figure 4: Same as in Fig. |^, but for the HO-improved predictions. 

Figure 5: LO predictions for the inclusive photoproduction of prompt Xcj mesons with 
radiative feed down to J/ip mesons in the form of (a) a integrated over 0.4 < z < 0.8 
and p\ > 1 GeV 2 as a function of W, (b) da/dp^ integrated over 0.4 < z < 0.9 and 
50 GeV < W < 180 GeV as a function of p\, and (c) da/dz integrated over p 2 T > 1 GeV 2 
and 50 GeV < W < 180 GeV as a function of z. The total results (solid lines) are 
built up by the direct-photon colour-octet (dashed lines), resolved-photon colour-singlet 
(dot-dashed lines), and resolved-photon colour-octet (dot-dot-dashed lines) contributions. 
Direct photoproduction is prohibited at LO in the CSM. 
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